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Abstract The common toad, Bu fo bu fo (Linnaeus, 1758), 
is widely distributed in Europe and parts of Africa. 
Previous studies of the genetic relationships among B. 
bufo populations have not included Asian populations. 
Here, we investigated the phylogenetic relationships 
of B. bufo, including a population from Xinjiang, 
China, using 722 bp of the mitochondrial cytochrome b 
sequence and identified six subclades and 85 haplotypes 
in populations from 22 countries. Phylogenetic analyses 
and pedigree differentiation revealed that the subclade 
containing the Xinjiang population has undergone 
a relatively recent expansion. Combining our results 
with those of previous studies, we found that the 
common toad population of Xinjiang, China, belongs 
to the European-Caucasian lineage and that its closest 
relatives based on phylogenetic relationships were 
eastern European populations. Our research contributes 
to expanding knowledge of the geographic distribution 
of B. bufo and illuminates the lineage and genetic 
relationships of the B. bufo population in Xinjiang, 
China. Future research should continue to update 
the geographic distribution of B. bufo and complete a 
genetic investigation of the full range of this species. 


Keywords distribution, genetic distance, genetic structure, 
haplotype, mitochondrial DNA, phylogenetic tree 


"Corresponding author: Prof. Yiming LI, from Institute of Zoology, 
Chinese Academy of Sciences, Beijing, China, with his research focusing 
on conservation biology and ecology of amphibians. 

E-mail: liym@ioz.ac.cn 

Received: 24 May 2020 Accepted: 11 September 2020 


1. Introduction 


The common toad, Bufo bufo (Linnaeus, 1758), is widespread 
in Europe, except Ireland, Iceland and some Mediterranean 
islands, and is also distributed in the western parts of North 
Asia and a small part of Northwest Africa (Lüscher et al., 
2001; Sillero et al., 2014). The evolutionary history of this 
species has been controversial, and its species delimitation and 
phylogeography has recently been revealed at the continental, 
national and regional scales (Arntzen et al., 2013; Cossu et al., 
2018; Garcia-Porta et al, 2012; Ozdemir et al, 2020; Recuero et al, 
2012; Tuncay et al, 2018). Garcia-Porta et al. (2012) studied 304 
samples from 164 populations covering most of the distribution 
range of the B. bufo species complex in the Western Palearctic. 
Studies show that the B. bufo species complex encompasses three 
main lineages: the Caspian lineage, the Iberian-African lineage, 
and the Europea n-Caucasian lineage. 

Although B. bu fo has been studied in areas where it is widely 
distributed, phylogeographic analyses of this species from some 
geographical areas, including Asia and specifically China, have 
never been included. B. bufo was discovered in Habahe County, 
Altay Prefecture, Xinjiang Uygur Autonomous Region, China, 
in late summer 2005 (Shi et al, 2005), and thus far, it is the only 
place in China where the common toad has been reported. 
However, there is limited information or historical data 
available to reveal the lineage and genetic relationships of the 
common toad in Xinjiang, China. 

B. bufo is a widespread and adaptable species present in 
coniferous, mixed and deciduous forests, groves, bushlands, 
meadows, arid areas, parks and gardens. It usually inhabits 
damp areas with dense vegetation, and large open areas are 
generally avoided. The species spawns and larval development 
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takes place in still waters and slow-moving parts of rivers 
and streams. It is present in many modified habitats (UCN, 
2020). In the Xinjiang region of China, basins are surrounded 
by high mountains, oasis and desert are intertwined, and the 
geographical environment of forests, grasslands and lakes 
(Wang et al, 2006) is suitable for the survival and reproduction 
of a common toad population. Habahe County, Altay, is at the 
northwesternmost edge of Xinjiang, bordering Russia in the 
north. It is located at the same latitude as France, Romania 
and other European countries, and the similar latitudes may 
cause these areas to enjoy similar climatic characteristics. This 
similar climate and suitable habitat are likely to have allowed 
the common toad to migrate from Europe to East Asia in the 
distant past. 

To investigate the phylogenetic information of the 
common toad found in Xinjiang, our laboratory working 
group conducted two scientific investigations in Altay 
Prefecture, Xinjiang, China, in 2014 and 2018, and successfully 
collected adult toad specimens. Previous researches has shown 
that genetic divergence can be used to determine whether 
individuals and populations belong to the same or separate 
evolutionary lineages (Ferguson, 2002) and can indicate the 
level of connectivity among populations (Blair et al., 2013; 
Dubey et al., 2011). Therefore, the informative mitochondrial 
cytochrome b (cytb) gene was selected to perform phylogenetic 
analysis and genetic structure analysis using new and 
previously published data. 


2. Materials and Methods 


2.1. Sampling Surveys were carried out in Xinjiang Uygur 
Autonomous Region in late summer of 2014 and 2018. 
Fifty-eight samples across three localities were used for the 
phylogenetic analysis. The third toe of the right hind foot from 
each postmetamorphosis common toad was clipped if they were 
available, and the tissue samples were preserved separately in 
95% ethanol and stored at -20°C in the laboratory. In this study, 
our laboratory extracted 58 common toad sequences from 
Xinjiang, China, and obtained all published sequences of B. bu fo 
among a total of 262 cytb sequences from NCBI, along with 
three outgroup sequences. Three sequences in the gene library 
had inconsistent lengths and were not included in the analysis. 
The locality and voucher data for sequences used in our 
analysis are given in Table Sl. The outgroups consisted of Bufo 
gargarizans gargarizans, Bufo bankorensis, and Bufo japonicus 
formosus, based on Recuero et al. (2012). 


2.2. DNA extraction and amplification In the laboratory, 
genomic DNA was extracted using the Universal Genomic 
DNA Kit (catalogue no. CW2298M; Beijing, CoWinBiotech 
Co., Ltd, Beijing, China) and following the manufacturer’s 
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instructions. To infer the putative geographic origin of the 
sampled individuals, we amplified one diagnostic mitochondrial 
marker, corresponding to a portion (722 bp) of the protein- 
encoding cytb region, using published primer pairs (Palumbi et 
al, 1991; Recuero et al, 2012). The amplification conditions were 
5 min at 94°C; 35 cycles of 30 s at 94°C, 30 s annealing at 48°C, 
and 30 s at 72°C; and a final 10 min at 72°C. The PCR products 
were then separated by electrophoresis on 2% agarose gels. The 
resulting PCR products were directly sequenced by Tsingke 
Company (Beijing, China) using the same primers used for 
amplification. The obtained sequences were compared with the 
available homologous sequences from GenBank (https://www. 
ncbinlmnih.gov/) using the Basic Local Alignment Search Tool 
(BLAST, http;//blast.ncbi.nIm.nih.gov/Blast.cgi) with the default 
parameters. Sequencing data are available at NCBI under 
accession numbers from М'Т499388 to MT499445 (Table S1). 


2.3. Data analysis The mitochondrial cytb gene sequences were 
aligned with reference sequences from Recuero and Arntzen 
using the ClustalW (Thompson et al, 1994) Multiple alignment 
program, which automatically applies gaps and sliding 
alignments, in BioEdit 7.2.5 (Hall, 1999). 

The haplotype distribution, number of all sequences, number 
of sequence sites in a sample (S), number of haplotypes (NOH), 
haplotype diversity (Hd) (Nei, 1987), sequence diversity (K) 
and nucleotide diversity (x) (Nei and Tajima, 1981) of each 
population were estimated using DnaSP ver. 6.13.03 (Rozas et al, 
2017). In addition, to study the degree of differentiation between 
populations and the amount of gene flow, DnaSP ver. 6.13.03 
was also used to calculate the differentiation parameter Fst 
between populations (Weir and Cockerham, 1984). Uncorrected 
pairwise sequence divergences (p-distances) between paired 
populations were calculated in MEGA v7 (Kumar et al, 2016). 

Arlequin 3.11 software was used to calculate the Tajima’s D 
(Tajima, 1989) and Fu's Fs (Fu and Li, 1993) of each population 
and perform a molecular variation analysis (AMOVA) among 
populations and clades to calculate the genetic structure and 
population genetic variation distribution (Excoffier and Lischer, 
2010). Arlequin was also used to calculate the raggedness index, 
R, (Ramos-Onsins and Rozas, 2002), Fu’s Fs value, Tajima’s D 
and sum of squared occurrences (SSD) value of each progressive 
branch to detect population expansion (Harpending, 1994). 
Based on the observed mismatch distribution and the expected 
value under rapid expansion, the SSD value showed a 
significant indication of population sta bility without expansion. 

Mismatch distribution analysis and estimation of 
parameters were performed using Arlequin 3.11 software 
(Excoffier and Lischer, 2010). The mismatch distribution graphs 
were drawn with DnaSP ver. 6.13.03 (Rozas et al. 2017). Self- 
testing of the saliency hypothesis was repeated 1000 times to 
infer dynamic changes in the historical population number. 
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By examining whether the mismatch distribution curve of 
the population was a single or multipeak type and whether 
it deviated from the neutral test, we could speculate whether 
the population had expanded in the past. If the general 
population has undergone expansion or sustained growth in 
the past, its mismatch distribution curve presents a single-peak 
Poisson distribution, and the Tajima’s D neutral test shows 
a significant deviation from neutral mutation (Slatkin and 
Hudson, 1991). However, when the population size has been 
stable, the mismatch distribution curve shows a multipeak 
distribution, and the Tajima’s D value is not significant 
(Schneider and Excoffier, 1999). Low Tajima’s D values and a 
neutral distribution of pairing differences indicate an ancient 
population expansion. 

The haplotypes of the cytb sequences were subjected to 
three types of phylogenetic analysis using MEGA v7 software 
(Kumar et al., 2016) and the Kimura-2-Parameter distance 
measurement (estimated gamma, «=0.32): the distance-based 
neighbor-joining (NJ) (Saitou and Nei, 1987), maximum 
likelihood (ML) (Kimura, 1980), and minimum evolution (ME) 
(Rzhetsky and Nei, 1992) methods. The robustness of each 
phylogeny was assessed by the percentage of 1000 bootstrap 
(BS) resamplings (Felsenstein, 1985). Phylogenetic analyses of the 
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haplotypes were also conducted using Bayesian inference (BI). 
BI analyses were performed in MrBayes 3.2.7 (Ronquist and 
Huelsenbeck, 2003). Each analysis used four heated Markov 
chains (using the default heating values) that were run for 1 
million generations. Trees were sampled every 1000 generations, 
and a consensus tree was calculated after omitting the first 1000 
trees as burn-in. 

Relationships among the haplotypes were visualized by 
reconstructing a median-joining network by Network 5.0.11 
(Bandelt et al, 1999). 


3. Results 


31. Haplotype diversity and genetic diversity Overall, a total 
of 320 sequences with a length of 722 bp were used (including 
three outgroup sequences) (Table Sl and Figure 1). For 317 
sequences of B. bufo, among the 722 loci, 625 were conserved 
and 97 were variant sites (including 58 simple information sites 
and 39 single base variant sites). 

The haplotype diversity (Hd), nucleotide diversity (x), and 
sequence diversity (K) of 22 different geographic populations 
are presented in Table 1 (outgroups excluded). Haplotype 
partitioning was performed for 317 cytb gene sequences, and a 
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Figure 1 B. bufo sampling localities included in the present study. Red triangles indicate the sampling localities in NCBI and other 


studies; blue triangles show the new found localities in China. 
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Table 1 Diversity indices and population expansion statistics of B. bufo. 
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total of 85 haplotype sequences were obtained, among which 
74 (87.06%) were unique and most appeared in only a single 
geographical population. A total of 12.94% of the haplotypes 
were shared and had a higher frequency. Of the 22 populations, 
17 contained haplotype polymorphisms. The Italian population 
contained 18 haplotypes followed by Turkey, while the 
Xinjiang population in China contained only two haplotypes, 
of which haplotype 2 was the unique haplotype, and haplotype 
1 contained the largest number of samples. 

Among all the populations, the sequence and nucleotide 
diversities were higher for the French population, Montenegro 
population, and Italian population than the other in situ 
populations, respectively. For China’s Xinjiang population, 
these two indicators were far lower than the other geographic 
populations at 0.189 and 0.00026. 


3.2. Genetic differentiation and genetic structure analyses 
Table 2 shows the statistical results of the differentiation index, 
Fst between 18 common toad populations; the four remaining 
populations provided no statistical results due to the small 
number of individuals. Overall, for the populations, the mean 
pairwise Fst varied from —0.139 to 0.914, and the analyses 
revealed significant genetic differentiation between most 
populations of B. bufo. The differentiation indices of Greece, 
Turkey and Italy compared to China were 0.804, 0.794, and 
0.713, respectively. These high Fst values indicated a long-term 
restriction of gene flow between these populations. China’s 
Xinjiang population showed the lowest differentiation index 
when compared with the population from Poland, and the 
indices for the populations from Croatia, Hungary, Slovakia, 
and Ukraine were also low. 

AMOVA was performed on all the clades and populations of 
the common toad (Tables ЗА—В). It can be seen from the results 
that the molecular variation was mainly between clades, which 
accounted for 82.37% of the total variation. The molecular 
variation within clades accounted for only 17.63% of the total 
variation. The differentiation index across all the clades was 
0.82375. After the permutation test, the P value was significant 
at < 0.01. A molecular variation analysis of the populations 
showed that the molecular variation was mainly between the 
populations, which accounted for 57.917% of the total variation. 
The molecular variation within the populations accounted for 
42.093% of the total variation. The differentiation index across 
the entire population was 0.57909. After the permutation test, 
the P value was significant at < 0.01. 


3.3. Phylogenetic analyses Analysis of the results of the four 
phylogenetic trees constructed with the toad haplotypes showed 
that the topological structures of the four phylogenetic trees 
and the bootstrap values of the clades were basically equivalent 
(Figure 2 and Figures SIA-C). A median-joining network 
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diagram of the toad haplotypes was constructed using network 
software with 11 intermediate vectors (Figure 3). The topology 
showed phylogenetic relationships that were basically consistent 
with those of the haplotype tree, with 85 haplotypes. 

According to the degree of clade support and the distribution 
of haplotypes, the common toad haplotype phylogenetic tree 
was from one same main clade divided into six subclades, Cladel 
to Clade6. These clades had obvious regional characteristics 
related to their distribution areas. 

Cladel was composed of 16 populations distributed in China 
and eastern Europe, and the Xinjiang population harbored a 
unique haplotype, Hap2. The 17 haplotypes in Clade2 were all 
distributed in Turkey, and one individual in Hapl7 was from 
Greece. Clade3 consisted of French and Italian populations, 
and their shared haplotype was Hap9. Clade5 mainly included 
the geographical population of Greece and its surrounding 
four countries. Clade4 and Clade6 both contained Italian 
populations. 


3.4. Pedigree Geography and Differentiation The results of 
the mismatch distribution analysis for each clade are shown in 
Figures 4A—F. Cladel, Clade2, Clade5and Clade6 each showed 
a single peak in the mismatch distribution curve, indicating 
that the populations in these areas have experienced expansion 
and rapid growth. The other branches showed bimodal or 
even multimodal distributions, indicating that the sizes of the 
populations distributed in these areas have remained relatively 
stable. 

The statistical results describing the differentiation and 
expansion of each clade of common toad are shown in Table 
4. The table shows that the SSD and R, tests did not reach 
significance in most clades (P > 0.05) except Clade2. These results 
are not completely consistent with the mismatch distribution 
map, which can only be explained if Cladel, Clade5 and Clade6 
have undergone rapid population expansion. Tajima's D and 
Fu's Fs based on the neutrality tests showed negative values 
for most clades at the clade level, and neither statistic reached 
significance (P > 005); these results could not confirm that the 
clade had undergone rapid expansion. In Cladel and Clade5, the 
value of Fu's Fs reached significance (P « 0.05), as did the value 
of Tajimas D (P « 0.05), indicating that the Cladel and Clade5 
populations expanded rapidly during their history which is 
consistent with the mismatch map. 

Tajima's D and Fu's Fs based on the neutrality tests were 
performed for each population, excluding the populations 
with fewer than five samples. The results showed that in 
most populations, these two indicators were negative but 
not significant (P » 0.05) (Table 1). In the Xinjiang population 
in China, neither Tajima's D value nor Fu's Fs indicated 
significance (P » 0.05); thus, this population may not have 


undergone rapid expansion. 
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Figure 2 Maximum Likelihood tree, representing relationship among haplotypes of B. bufo. 
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Table 3A Results of analysis of molecular variance (AMOVA) between clades of B. bufo. 


Source of variation d.f. Sum of squares Variance components Percentage variation 
Among clade 5 965.474 4.48097 Va 827 
Within clade 311 298.176 0.95877 Vb 17.63 

Total 316 1263.65 5.43973 


Fixation index 


Fst: 0.82375 


Table 3B Results of analysis of molecular variance (AMOVA) between populations of B. bufo. 


Source of variation d.f. Sum of squares Variance components Percentage variation 
Among population 21 726.889 2.50328 Va Al 
Within population 295 536.761 1.81953 Vb 42.09 
Total 316 1263.65 4.32281 
Fixation index Fst: 0.57909 
Table 4 Results of the statistical tests of neutrality and mismatch distribution of subclades of B. bufo. 
T 0 
Clade N NOH SSD R, Tamjia’s D Fu’s Fs 
(С1=95%) (С1=95%) 
0.5 1.06831 0.000932 0.02665 -2.37134 —28.37072 
Cladel 145 38 
(0.46709-5.41881) (0.00072-2.08568) (P=0.762) (P=0.855) (P<0.001) (P<0.001) 
3.3 0.00072 0.02451 0.07602 0.90863 —5.41644 
Clade2 66 17 
(1.76309-4.34540) (0.00072-1.59289) (P=0.012) (P=0.025) (P=0.180) (P=0.022) 
0.5 255563) 0.04837 0.08721 0.09233 0.10908 
Clade3 19 7 
(0.85484-8.90724) (0.00072—4.15897) (P=0.110) (P=0.326) (P=0.584) (P=0.557) 
2 0.00072 0.01939 0.26915 -1.25355 -1.13718 
Clade4 31 6 
(0.00001-4.76180) (0.00072-0.35158) 0.608 (P=0.753) (P=0.121) (P=0.215) 
1.4 0.01 0.00217 0.06129 -1.98343 —3.94238 
Clade5 12 7 
(0.87200-2.80740) (0.00072-1.08754) (P=0.636) (P=0.456) (P=0.006) (P=0.002) 
2:7 0.00072 0.00295 0.02161 –0.88266 -2.80691 
Clade6 28 10 
(1.41420-4.41574) | (0.00072- 1.44821) (P=0.809) (P=0.924) (P=0.215) (P=0.066) 


N: number of sequences; NOH: number of haplotypes; т: time in number of generations elapsed since the rapid expansion episode; 0: mutation 
parameter; SSD: sum of squared deviations; R2: raggedness index; CI: confidence interval. 


4. Discussion 


Our study reports the first use of mitochondrial cytb data 
to explore the lineage and genetic structure of B. bufo from 
Xinjiang, China, in comparison with other samples from 21 
different geographical areas. Based on the haplotype network 
and haplotype evolution tree, the common toad populations 
from the 22 different regions belong to the same lineage, 
consistent with the results of Garcia-Porta et al. (2012). China’s 
Xinjiang population, located in Cladel, emerged late in the 
differentiation process and belongs to the European-Caucasian 
lineage (Garcia-Porta et al., 2012). Based on the six subclades, 
the Xinjiang population of China has a close phylogenetic 
relationship with some Eastern European populations. 

The level of genetic diversity is the product of long-term 


evolution and an important parameter to evaluate the present 
situation of biological resources. The higher the genetic diversity, 
that is, the richer the genetic variation, the stronger is the 
adaptability to the environment and survival competitiveness. 
The nucleotide diversity index and haplotype diversity index 
are important indexes to measure the genetic diversity of a 
species population. Grant and Bowen proposed that when 
a population has a high h and a low zx (for example, h > 05, 
п < 0.005 or 0.5%), it can be speculated that the population 
may have been subjected to a bottleneck effect followed by 
rapid population growth in the past. Haplotype diversity and 
nucleotide diversity analysis (Grant and Bowen, 1998) showed 
that the populations of Hungary, Romania, Russia and other 
eastern European regions may have experienced bottlenecks 
followed by historical events of rapid population growth. The 
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Figure 4 Mismatch distributions within Cladel-Clade6 (A-F) of B. bufo. Dotted lines show the observed distribution of mismatches, 
and solid lines show the expected distribution under an expansion model. 


AMOVA and statistical results for the mismatch distribution 
of all the subclades of B. bufo also support the inference that 
the populations in Cladel underwent expansion and rapid 
growth in the past. Cladel included populations distributed in 
most countries in eastern Europe, a small part of the Russian 
population and the entire Chinese Xinjiang population. The 
Fst results revealed the existing genetic differentiation of this 
species. A lack of gene flow was observed between the common 
toad population of Xinjiang and the populations from Italy, 
Turkey, and Greece. In contrast, the genetic differentia tion 


difference was smaller between the Xinjiang population and 


eastern Europe population. 

Based on the above data, it is speculated that the populations 
from Eastern European countries likely rapidly expanded 
rapidly in historical time. General speaking, the high level 
of genetic differentiation between populations was due to a 
lack or low level of gene flow, especially between populations 
located far away in space or separated by natural barriers 
(Weerachai et al., 2019). For example, due to the barrier of 
the Alps, the Italian population in this study demonstrated a 
very weak degree of gene communication with the external 


population, thus forming a significant evolutionary branch in 
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the phylogenetic relationship and demonstrating great genetic 
differentiation from the other populations. The geographical 
area in which the population from the Eastern European 
countries is located belongs to the plain area of Eastern Europe, 
which has a relatively low average elevation and no high 
mountains and rivers, features that are conducive to the spread 
of amphibian species. 

In addition, climatic fluctuations and periodic glacial action 
in the Pleistocene contributed to the current patterns of genetic 
variation in biological species in Europe and North America 
(Hewitt et al, 1996; Hewitt et al, 2000). During the glacial period, 
especially during the last ice age (LGM), the climate was dry, 
the temperature was low, and large-scale glaciers were present 
globally (Gasse et al., 2000; Clark et al., 2009). Quaternary 
glaciers have a significant impact on the vast areas of Europe, 
and the southern edge of the European ice cap can reach 50" N. 
Many species experience severe habitat shrinkage and continue 
to survive only in restricted geographical shelters. During warm 
interglacial periods, shrinking ice caps create new suitable 
habitats, and species begin to expand from shelters. 

It can be inferred that a common toad population or a 
few individuals probably experienced bottleneck effects due 
to extreme climatic conditions. Subsequently, the population 
underwent a rapid expansion process, spreading eastward 
from Eastern Europe, and finally entering Xinjiang from 
the northern border. However, further research is needed to 
support this inference. 

However, according to the national amphibian survey, the 
common toad is not distributed in any part of China outside 
Xinjiang Uygur Autonomous Region (Fei et al., 2012; Shi et al., 
2005). Neutrality tests based on population Tajima's D and Fu's 
Fs showed that these two indicators were not significant (P » 
0.05) in the Xinjiang population, indicating that this population 
may not have undergone a rapid expansion. This lack of 
expansion may have been caused by the unique diversity 
and particular traits of the complex Xinjiang ecosystem; the 
presence of the Gobi desert and other arid regions could have 
prevented the common toad population from spreading further 
(Wang et al., 2006). The haplotype network graph shows that 
the Xinjiang sample in China has a unique haplotype (Hap2), 
indicating that this population has survived and evolved for a 
long time in Xinjiang, China, and is isolated from its original 
parent population, resulting in mutation and genetic drift. 

To understand the origin and future diffusion patterns 
of a species, it is very important to investigate and study the 
distribution regions and phylogenetic relationships within 
the species. Although the mitochondrial partial sequence data 
used in this study may not be sufficient to accurately reflect 
the whole evolutionary history of B. bufo, it can at least 
reveal the general framework of the evolutionary history 
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of this species. Therefore, in the future, we will exam the 
complete mitochondrial genome sequence to more reveal the 
evolutionary history of species. Further studies of geography, 
genetics and molecular biology should also be conducted to 
provide effective conservation for more species. 
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Figure S1 Phylogenetic tree, representing relationship among haplotypes of B. bufo. A: Minimum Evolution tree; B: Neighbor-Joining 
tree; C: Bayesian tree. 


Table 51 Sequence information of B. bufo from 22 countries in this study. The ACCN numbers marked in red are not included in the analysis 
because they are too short. The ACCN numbers marked in blue are three outgroup sequences. JN from Recuero et al. (2012), MF from Arntzen 
et al. (2017), MG from Cossu et al. (2018), MN from Ozdemir et al. (2020), KX from Tuncay et al. (2018) and MT from this research. 


ACCN Latitude (N) Longitude CE) Country 


JN647248 52.30426 5.825339 Netherlands 


JN647250 44.60739 20.57056 Serbia 


JN647252 44.68219 20.56553 Serbia 


JN647254 44.68219 20.56553 Serbia 


JN647256 43.77905 19.99715 Serbia 


JN647258 43.77905 19.99715 Serbia 


N647260 43.78008 19.98942 Serbia 


= 


N647262 39.72456 21.75401 Greece 


= 


N647264 43.4219 20.37977 Serbia 


= 


N647266 42.89633 19.64393 Montenegro 


= 


N647268 47.00947 17.63034 Hungary 


J 


N647270 46.09315 18.14138 Hungary 


— 


N647272 42.61611 21.43511 Serbia 


== 


JN647274 42.43333 21.90472 Serbia 


JN647276 42.45083 21.88917 Serbia 


JN647278 43.77905 19.99715 Serbia 


JN647280 64.1 18.3833 Sweden 


JN647282 64.1 18.3833 Sweden 


JN647284 64.1 18.3833 Sweden 


JN647286 55.71089 36.77122 Russia 


JN647288 56.84039 60.55564 Russia 


JN647290 41.18069 35.77406 Turkey 


JN647292 56.53 31.53514 Russia 


JN647294 53.34931 35.55847 Russia 


JN647296 66.53567 33.16139 Russia 


JN647298 55.86925 49.14461 Russia 


JN647300 45.23331 28.30869 Romania 


Continued Table 81 


ACCN Latitude CN) Longitude (E) Country 


JN647302 39.56042 21.37186 Greece 


JN647304 41.38959 24.63762 Greece 


JN647306 38.56438 22.34908 Greece 


JN647308 47.8669 24.5154 Romania 


JN647310 46.7827 25.7838 Romania 


JN647312 45.4424 24.6061 Romania 


JN647314 49.535 3.677778 France 


JN647316 50.01667 3.875 France 


JN647318 49.5169 3.646967 France 


JN647320 45.63102 16.44545 Croatia 


N647322 49.6261 36.28553 Ukraine 


N647324 49.61971 36.3139 Ukraine 


N647326 50.82114 1.602306 France 


N647330 57.73333 11.91667 Sweden 


N647332 57.73333 11.91667 Sweden 


= 


N647334 50.69778 7.128333 Germany 


N647336 51.7891 15.7282 Poland 


— 


N647338 50.20013 20.79382 Poland 


= 


N647340 50.26361 5.368611 Belgium 


= 


N647342 42.36956 18.89072 Montenegro 


— 


N647344 47.35313 17.77907 Hungary 


JN647346 47.26872 17.6948 Hungary 


JN647348 46.96708 16.50817 Hungary 


JN647350 46.9092 17.84924 Hungary 


JN647352 46.21503 18.61146 Hungary 


JN647354 45.23903 28.297 Romania 


JN647356 45.23903 28.297 Romania 


JN647358 41.39108 32.81772 Turkey 


ACCN Latitude CN) Longitude (E) 


JN647360 40.82817 36.60161 


JN647362 38.73956 16.23598 


JN647364 39.16737 16.07913 


JN647366 42.2699 12.06063 


12.761 


JN647368 41.7501 


JN647370 41.7501 12.761 


JN647372 39.566 16.05428 


JN647374 39.25795 16.09573 


JN647376 39.25795 16.09573 


JN647378 41.47082 12.71834 


15.84592 


JN647380 38.181 


JN647382 44.2483 11.34623 


JN647384 41.88985 13.17065 


JN647386 41.88985 13.17065 


JN647388 39.94193 15.80598 


JN647390 39.94193 15.80598 


JN647392 41.26808 13.04561 


13.04561 


JN647394 41.26808 


JN647396 42.14047 12.09707 


JN647398 42.14047 12.09707 


JN647400 39.80038 15.90794 Italy 


JN647402 36.92889 14.67362 Italy 


JN647404 36.92889 14.67362 
JN647406 38.47813 16.46922 


JN647408 41.84114 13.03899 


JN647410 44.03867 10.86172 


JN647412 42.27485 12.92563 


JN647414 42.27485 12.92563 Italy 


Continued Table $1 


ACCN Latitude CN) Longitude (E) Country 


JN647416 41.01962 15.00642 Italy 


JN647418 41.01962 15.00642 Italy 


MF461245 46.577 4.016 France 


MF461247 46.577 4.016 France 


MF461252 46.273 6.771 France 


MF461254 43.119 6.318 France 


MF461256 43.119 6.318 France 


MF461261 46.358 5.556 France 


MF461263 46.358 5.556 France 


MF461265 43.891 7.486 France 


MF461267 45.13 5.588 France 


MF461269 45.13 5.588 France 


MF461271 45.219 7.541 Italy 


MG199042 40.53567 9.08895 Italy 


JN653295 25.02069 121.5601 


5 


MT499388 48.06006 86.39969 


MT499390 48.06006 86.39969 


MT499392 48.06006 86.39969 


MT499394 48.06006 86.39969 


MT499396 48.06006 86.39969 China 


MT499398 48.06075 86.38649 China 


MT499400 48.06075 86.38649 


MT499402 48.06075 86.38649 


MT499404 48.06075 86.38649 


MT499406 48.06075 86.38649 China 


MT499408 48.06075 86.38649 China 


MT499410 48.06075 86.38649 China 


MT499412 47.99625 86.28246 China 


Continued Table S1 


ACCN Latitude (N) Longitude (E) Country 


MT499414 47.99625 86.28246 China 


MT499416 47.99625 86.28246 China 


MT499418 47.99625 86.28246 China 


MT499420 47.99625 86.28246 China 


MT499422 47.99625 86.28246 China 


MT499424 47.99625 86.28246 China 


MT499426 47.99625 86.28246 China 


MT499428 47.99625 86.28246 China 


MT499430 47.42656 85.70167 China 


MT499432 47.42958 86.50242 China 


MT499434 47.42958 86.50242 China 


MT499436 48.0314 86.48003 China 


MT499438 48.0314 86.48003 China 


MT499440 48.0314 86.48003 China 


MT499442 48.0314 86.48003 China 


MT499444 48.0314 86.48003 China 


MN432924 37.459857 30.904502 Turkey 


MN432926 37.459857 30.904502 Turkey 


MN432922 37.766994 28.972456 Turkey 


MN432978 37.936607 27.422264 Turkey 


MN432980 37.936607 27.422264 Turkey 


MN432948 38.506016 31.311516 Turkey 


MN432981 39.864236 29.6401423 Turkey 


MN432983 39.864236 29.6401423 Turkey 


MN432970 40.055243 26.4214873 Turkey 


MN432975 40.610907 29.070776 Turkey 


MN432977 40.6 10907 29.070776 Turkey 


MN432960 41.929597 27.366404 Turkey 


Continued Table 81 


ACCN Latitude (N) Longitude (E) Country 


MN432984 41.22 28.96 Turkey 


MN432986 41.22 28.96 Turkey 


MN432954 40.568289 30.004928 Turkey 


MN432918 40.606408 31.2866436 Turkey 


MN432920 40.606408 31.2866436 Turkey 


MN432956 41.430675 32.557278 Turkey 


MN432972 41.426331 33.008059 Turkey 


MN432974 41.426331 33.008059 Turkey 


MN432958 41.402751 33.697383 Turkey 


MN432965 42.0245258 34.9196832 Turkey 


MN432967 42.0245258 34.9196832 Turkey 


MN432951 41.669594 36.035401 Turkey 


MN432938 40.75732 37.429212 Turkey 


MN432940 40.75732 37.429212 Turkey 


= 


N432935 40.648342 38.45252 Turkey 


= 


N432941 40.617968 40.29556 Turkey 


Es 


N432943 40.617968 40.29556 Turkey 


Es 


N432932 41.027744 41.003802 Turkey 


KX230485 60.1656 5.9369 Norway 


KX230487 60.0172 5.4853 Norway 


KX230489 58.0692 7.9939 Norway 


63.4841 8.8846 Norway 


KX230483 58.0692 7.9939 Norway 


58.9488 5.8904 Norway 


59.0175 5.7803 Norway 


60.3035 5.5061 Norway 


61.8039 4.954 Norway 


63.5165 9.0323 Norway 


Continued Table 81 
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Latitude CN) Longitude (E) Country 


63.4834 8.7855 Norway 


~ 
o 
EN 
со 
~ 


58.0692 7.9939 Norway 


58.7855 5.7348 Norway 


59.7037 5.3937 Norway 


59.7772 5.4685 Norway 


59.7841 5.2553 Norway 


60.1656 5.9369 Norway 
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